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The Green’s-function technique is utilized in the determination of the field distribution of an infinitely 
long current carrying conductor of arbitrary but constant cross section above a superconducting ground 
plane of finite thickness. It is assumed that the superconductor can be described by the phenomenological 
London equations. The integral expressions that are obtained are solved analytically for a few special cases 
of interest. Under conditions that are often encountered in a physical system, a modified image method can 
be utilized in order to calculate the field distribution to within 2% of the computer solution, 


INTRODUCTION 

D UE to the rapid advance of the technology of super- 
conductive devices for computer applications, a 
great deal of interest has been expressed in developing 
high-speed switching components. 1 A device whose 
basic structure is that of a thin superconducting film 
above a superconducting ground plane exhibits the 
desired switching speed characteristics. 2 * 3 The basis of 
any electromagnetic analysis of this type has been the 
assumption that the return current in the ground plane 
is localized under the current carrying film, effectively 
forming a strip transmission line whose characteristics 
are then calculated. 4,5 The assumption (henceforth 
called the strip-line assumption) is made for materials 
whose film width is much larger than the separation 
between the film and ground plane. 

It is often necessary to determine the field distribu- 
tion for materials that do not fit the requirement of the 
strip line assumption. In this paper, the Greeks func- 
tion technique is used in deriving a general expression 
for the field of a conductor that is carrying current 
above a superconducting ground plane. The strip line 
assumption is verified and a few cases of particular 
interest are examined. A modified image technique is 
developed that is useful for determining the field of the 
current carrying conductor if some geometrical condi- 
tions that are often encountered in practice are satisfied. 

CHARACTERIZATION OF THE SUPERCONDUCTING 
GROUND PLANE 


where /J is the reciprocal London-penetration depth and 
B = /x 0 H. (The rationalized mks system of units is used 
throughout this paper.) It can be shown from Eqs. 

(1) and (2) that in the Coulomb gauge under static 
conditions 

Vx(A+moT 2 J)=0, (3a) 

T-(A+moT 2 J)=0, (3b) 

and 

V 2 J=0 2 J. (4) 

The following statements are considered to be valid 
in the following analysis (Fig. 1) : (1) J in the ground 
plane has only one component (in the z direction); 

(2) J in the ground plane is bounded in space (this is 
justified later). It is assumed that the current density 
in S' is known. In a suitable gauge 

A=-/xod- 2 J (5) 

and 

V 2 A=0 2 A. (6) 

A straightforward argument can be used to prove 

that J is uniquely determined for a given total ground 

plane current. 

STATEMENT OF PROBLEM 

Consider the situation in which a current carrying 
source of arbitrary cross section is placed above a super- 
conducting ground plane (Fig. 1). A fraction £=/i// 0 
of the current is returned through the ground plane and 
(1 — £) of the current is returned through the wire above 


It is assumed that the properties of the ground plane 
can be described by the London equations 


E=Mo(/ 8) 2 (dJ /dt), 

(1) 

H = — (/})-* vxj, 

(2) 
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the conductor. It is assumed that there is no geometric 
variation in the z direction so that d/dz of the field 
quantities are zero. 

The source current density J s(x\y'), which has only 
a z component, can be separated into a symmetric and 
antisymmetric part (with respect to y') called J* and J a , 
respectively. 

J'&,rt=hlh&J)+W, -/)], (7a) 

= W, -?')], (7b) 

where 


L 


j a (*'yMS'=/ 0 . 


( 8 ) 


The Green's function diagram, which corresponds to 
Fig. 1, is shown in Fig. 2. Due to the separation defined 
by Eqs. (7a) and (7b) symmetric and antisymmetric 
pairs of delta function sources are used to generate the 
Green's function of the system. For simplicity consider 
the case for £ = 1 first. Since free space can be described 
by the equation 

V*A«-/ioJ, (9) 

the solutions to Eqs. (5) and (9) in the regions shown 
in Fig. 2 are as follows: (A Z8 (t/t') and A 2a (r/T r ) are the 
symmetric and antisymmetric (Green's) vector poten- 
tials, respectively). 


Region 7, 2: 

’A za ( r/r')i .2 
.(i/Oi* 


pM r/r')i. 2 l_ rrC u (k)coskyn 

L.4 20 (r/r')i.2-l Jo Lcia(fc) sinfryJ 


e~ k *dk\ (10a) 


Region 3: 

r4„(r/r') 3 ~l r*rcosky~] 
L.4 ra (r/r') 3 J Jo LsinfcyJ 



Region 4: 

rA z ,{T/t')n r-rcos^-1 
LA 2 „(r/r')iJ Jo L sin&yJ 


(10b) 


(10c) 


Region 5: 


[7 c '‘ C0Sky V dk 

L„(r/r')J Jo\C, a sinky/ 
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Fig. 2. Green’s function diagram. 


The solution e kx was rejected for x>x\ and the solu- 
tion e~ kx was rejected for #< — </, since it is assumed that 
the current density in the ground plane is bounded (in 
the x,y plane) and that the Green's vector potential is 
zero at x,y=<*>. This is reasonable since the ground 
plane in any physical situation is not infinitely wide 
but can indeed be very wide in comparison to the 
source conductor dimension and x\ A similar statement 
can be made regarding boundedness in the z direction. 
Therefore, the boundedness statement given in the last 
section is justified by the above physical argument. 

In order to evaluate the C's, it is necessary to realize 
that the x component of B is continuous at x=x\ 0, 
— d, while the y component of H is continuous at # = 0, 
— d and is discontinuous at x=x f . The five equations 
describing the continuity of field are simply algebraic in 
form and are not explicitly shown here. The discon- 
tinuity of H y at x=x f is not quite so obvious and is 
derived here. At x=x' 

H y ( region 1, 2) — /^(region 3) 

~6(y-y')±5(y+y'). (11) 


The signs on the right refer to the symmetric and anti- 
symmetric cases, respectively. 5(yzhy') is the Dirac 
delta function. Eq. (11) becomes 



/cos ky\ 

(k)-C 2 (k)^ke- kx '+C z (k)ke- kx '}[ m )dk 

\ sin ky) 


=mo C«(y-y)±8(y+/)]. (12) 


Using the Fourier integral theorem yields 

2ju 0 /cos ky\ 

[C,(*)-C,(*)>-^+C, (*)*•'-— ( J. (13) 

kir \sin ky / 


Simultaneous solution of Eq. (13) and the five con- 
tinuity equations yields the values for the C's. If 5 is 
defined as cos ky' in the symmetric case and sinky' in 
the antisymmetric case and 

a ± =k±(l3*+k 2 y) 

y± = e ±(eh-k,)U L (i 3a ) 

5±= lie 2 **' 
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then 

5 2/io r a+—k 

Ci= e t(x '~ d> \ (S-) (Y + aq_+7_a_) 

A ir L k 


+ (5 + )(7_a_— 7+“+) 

-5 2/zo 

C 2 = e *(x'-.)[ (a _)( a+ )(7 + _ 7 _)] 

A kir 

Cz= (Sno/kT)e~ kx ' 

Q= ± (4Smo/ A Tr)e k(x '- d) a±y ± 

C 6 = (— S/A)(8no/Tr)e kx '(a+— k), 


(13b) 


where 

A = + 2e* (2x '“ d) [(Y+) (<*+) 2 - (7-) (a-) 2 ]. 

Once the Green's potentials have been found, it is a 
simple matter to express the vector potential for the 
configuration shown in Fig. 1 (remembering that the 
case for £= 1 is under consideration) : 

Regions I, II, IV, V, VI: 


A z (x,y) = j J [_J,(x',y')A„(r/r') 

iS’ 

+J a {x'y')A 2a (r/r')']dx'dy', (14) 


(\S r refers to an integration over the source area defined 
by />0) where the primed integration is over the 
source cross section. A 28 (t/t') and Aza(t/i') are the 
appropriate vector potentials for the region under 
consideration. 


Region III: 

A,(x,y) = J j lJ,(x',y')A !t ( r/r')i 

+J a {x',y')A J<x (r/r')i ~\dx'dy' 


Si 


+11 U'W,y' 


)/l„(r/r')3 


s 2 


~VJa (%',y')A za (r'/ r') z^dx'dy', (15) 
vvhere S^ is the portion of the source region in which 
x<x\ 
v'>0, 

and S? is that portion in which 
x>x', 
y'> 0. 


(16a) 


(16b) 


In the above equations, r' and r are the source and field 
points, respectively. 

It can now be shown that the total current in the 
ground plane is — 7 0 . From Eqs. (5) and (14) it is seen 
that the current in the ground plane is 

'--=////* * (x V )A 28 (r/x')i~]dx'dy'dxdy, (17) 

S G P *S' 


where the unprimed coordinates refer to a point in the 
ground plane. (It can be shown by direct integration 
that the antisymmetric current density does not con- 
tribute to the net ground-plane current.) Integrating 
over y first yields 


/= 



«(*)/.(* V) 


X[C 4s e^^^ x +C 5 ^- ( ^ Ar2) ^]^x'(/y'7a:. (18) 


Integrating ovep k and x, respectively, yields 


r - 2 // 


J 8 (x'y')dx'dy'=-I 0 , 


*S' 


(19) 


which completes the proof. 

PARTIAL CURRENT RETURN IN THE 
GROUND PLANE (0<£<1) 

This case can be treated as the superposition of two 
situations. In the first situation, all the source current 
/o is returned through the ground plane. In the second 
situation there is no source conductor, and the return 
wire carries current (/i— 7 0 ), which is returned through 
the ground plane. 

Using this superposition the Green's vector potential 
can be found directly by considering the case for £= 1. 

Define constants D(k) that are related to C(k) by 
the relations 

D 8 (k) = -(l-Z)C 8 (k)/2 (20a) 





(and replace x' in the C 8 (k) expressions by x") 

D a (k) = 0. (20b) 


Thus, the Green's potentials are: 

Region 1: 

A 2 (r/r') l = A 2 (r/r') l \^ l + I Di(k) coskye~ kx dk ; (21a) 
Jo 

Region 2: 

i4,(r/r') 2 =^4*(r/r')ih-i 

+ / cosky\_D 2 {k)e~ kx + Dz{k)e kx ~]dk\ (21b) 
Jo 
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Region 3: 

^ z (r/r')3=.4 £ (r/r') 8 |{-i 

+ f cosky\_Di{k)e~ kz +D%{k)e kx ~\dk] (21c) 
Jo 

Region 4: 

.4 J (r/r') 4 =.4 2 (r/r) 4 |{_i+ 1 cosky 

Jo 

XlDiWe^+W'+DsWe-W+W'^dk; (21d) 

Region 5: 

i4»(r/r , )s=i4 l (r/r , )s|{_i+ f D t (k) coskye kx dk. (21e) 
Jo 

For many cases of practical interest, the return wire 
is far from the ground plane and the source conductor. 
It can be verified that 

lim D.(*) = -C(0)(l-€)[*(*)/«(0)] 

x"— +oo 


can be used to advantage in calculating the inductance 
of the structure. 

^ The Green’s vector potential is now completely 
specified for arbitrary £ so that it is now appropriate to 
examine a few source conductor geometries of general 
interest (it is assumed that the return wire is essentially 
at infinity for the remainder of the paper). 

SPECIAL CASES OF INTEREST 
Case 1 

h(x\y r ) = J L d(x f -a). 

This source distribution violates the requirement that 
/ o be finite so that Eq. (22b) cannot be applied here. 
Although the fields from this distribution can be more 
easily obtained by a one-dimensional analysis, they can 
also be obtained from the above equations to demon- 
strate the validity of the analysis. Using Eq. (14) with 
Eqs. [21 (a-e)] for the Green’s vector potential and 
performing the curl operation yield the following 
relationships: 


n=2, •••, 6. (22a) 

Therefore, 

lim £ 4 t (r/r , ) B |0<f<l = ^*(r/r')l{-i»=2,- • ,6. (22b) 


The net current in the ground plane can be found in 
a manner identical to Eqs. (17) to (19) and can be 
shown to be —I i. This means that although the ground 
plane only returns a fraction of the source current, the 
vector potential is identical to the case for £=1 if the 
return wire is far from the source and the ground plane. 
If one remembers the superposition that is the basis of 
Eqs. (21a)-(e), it is apparent that if the ground plane 
carries a current —I i, — 7 0 of the current follows the 
distribution given by Eq. (17), while (7 0 — 7i) is uni- 
formly distributed over the entire ground plane cross 
section. The current density of the uniform distribution 
is infinitely small so that it does not contribute to the 
vector potential, however, its integral over the cross 
section of the ground plane is still (7 0 — /i). It should 
also be noted that Eq. (22a) is valid only if / 0 is finite. 
A geometry that violates this restriction is examined 
below. 

The above discussion is not as remote as it may seem. 
Consider a ground plane of finite width W Q . Assume 
that the source conductor has dimensions small com- 
pared to W 0 and is close to the ground plane surface. 
Also assume that the return wire is far from the source 
conductor and ground plane. However, stipulate that 
the perpendicular distance from the return wire to the 
ground plane is also much less than W g . For this case, 
Eq. (22b) is satisfied in the neighborhood of the source 
conductor. For certain geometries (discussed later) the 
fields in the neighborhood of the source are very large 
in comparison to the fields elsewhere so that Eq. (22b) 


x>a 
0 <x<a 

—d<x < 0 


Hy=J L { 1-0/2, (23a) 

H v =-J l ( l+Q/2, (23b) 

M 1+g) [ r(l-g)- (1+0 coshftq 

2 IL (l+£) sinh/3rf J 


X sinhjS.t— cosh/Sx } , 

x<-d H y =-J L { l-£)/2. 

Case 2 

\J L b(x’-a) \y’\<W 
0 \yf\>W 


Js(x',y') = 


(23c) 

(23d) 


(24) 


The current per unit thickness in the ground plane, 
which is contained in a width 21, is 

£2 -Z ~x /.TF -x 

Ji(x)= / / / / a) cosky 

Mo J-i J o J o J -x 


X[C 4 ,(ft)^* 8)l *+C5,(fc)^ (25) 


Integration with respect to x', y', and y, respectively, 
yields 


/*(*)=- 


2/3 2 r- r 

— / 4 

no J o l 


Ci 8 (k,x' = cl) 


. e {p+&)\x 


cos ky' 

* I 

cos ky' J k 


C ba (k,x' = a) “| sin kl sinHF 

-€-<*+»** I dk. (26) 


Introducing the change of variables k' = kd yields for 
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the case PF»/ 3 ~ 1 , a 

7i(*)« 

rfTT L e+w-e-o* J 

r x sin(k'l/d) s\nk’ (W/d) 

X dk'. (27) 

J o *'(//<*) *'<wy<o 

The current density is given by 

dJi(x)/d( 21) 

so that for KW 


r-^^d+z) g— 0(d+X)-| 


J,(x)~-I 3J l \ . 

(28) 

L e+t d -e-t d J 


When l > W 



(29) 


which verifies the strip-line assumption discussed in the 
introduction. 



Fig. 4. Ground-plane 
current density. 


MODIFIED IMAGE TECHNIQUE 

There is a large class of problems that can be solved 
without resorting to manipulations such as were shown 
in case 1 and 2 by using a modified image technique 
that is now discussed. For simplicity, consider the 
Green’s-vector potential that is generated by a single 
delta-function source of unit strength inside the source 
conductor at x' = a. There is no loss in generality by 
taking y' = 0. 

If di (see Fig. 1) is much larger than p~ l then a/p~ l y>l, 
since a is the x coordinate of an arbitrary delta-function 
source inside the source conductor. If a/P~ l ^>l the 
Green’s current density in the ground plane is [Eqs. (5) 
and (14) ; a factor of ^ must be introduced into Eq. (14), 
since here a single delta-function source is being con- 
sidered (rather than a symmetrical pair)] 

1 P cosh {[(#/</) + 1 ~]pd) 

J *(r/r') « ; 

7r d sinh pd 

r x ky 

X / cos — e~ ka,d dk y (30) 
Jo d 

which reduces to 


1 a cosh[(#/d) + ll&/ 

r/r') « p . 

7 r a 2 +y 2 sinh pd 


(31) 


The factor — a/7r(a 2 +y 2 ) is the current per unit width 
at the #=0 surface of the ground plane which is found 
by replacing the ground plane by the image of the delta 
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function source. Since Eq. (31) is valid for every delta- 
function source that makes up the source conductor, the 
following general rule can be established. Given a source 
conductor of arbitrary cross section with the specifica- 
tion that the current density in the ground plane 


can be found by first finding the current per unit width 
0) by using the method of images and then 
multiplying /u,(x=0) by cosh{[(#/d)+l]/3d}/sinh fid. 

Equation (31) reduces to the perfect conductivity 
limit (/3 —1 — > 0) as can be seen by considering the 
following : 


lim— (I/7 r)fieP x [_a/ (& 2 +y 2 )]. 

wo 

(32a) 

It is interesting to note that 


r 00 x—Q 
lim \ 

l 0 x<0, 

and 

fO 

(32b) 

lim / fieP x dx= 1. 

J-d 

(32c) 

Thus, for #<0 


lim fieP x =5(x) 

/S— >00 

(33) 

and 


r/r') = - (l/7r)[a/(a 2 +;y 2 )]S(*), 

(34) 


which is the correct result if the ground plane is a 
perfect conductor. 

It is now instructive to examine the numerical solu- 
tion to the Green’s current density in the ground plane 
for the case of a unit delta-function source at x' = a and 
y f — 0. The_ dimensionless expression for the current 
density [from Eqs. (5) and (14)] is 


where 


r 30 (Pd) 2 /y \ 

J z (x,y)d 2 = - / iifc— — cos(~k ){A ie u*M+ Ut(^) 2 +^i + ^ 2e -[(^)+m^)2 + / c 2 ]i}j 

JO 7T \d / 

\A J ~ ^[owh-a*] »{fe+|:0^ 2 4-* 2 ]*}2— 


(35a) 


(35b) 


Although the left side of Eq. (35a) does not appear 
to be dimensionless, it should be remembered that a 
unit strength source is under consideration. Equation 
(35) was solved by a modified Simpson’s rule and an 
IBM 7090 computer. It was found that the integrand 
converged rapidly to zero so that infinity could be 
replaced by a finite number N (the highest N that was 
used was 81). The results are given in Figs. 3 to 6. It 
should be noted by comparing Eq. (31) with Figs. 3 and 
5 that the modified image technique developed in this 
paper is valid to within 2%. The numerical results for 


the Green’s current density can be used to evaluate 
the current density in the ground plane for an arbitrary 
source conductor by using approximation techniques 
for cases in which the modified image technique is 
not valid. 
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